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Abstract. We study the correlations with inclination of 
id-band disk and bulge structural parameters and I — H 
colour profiles for a sample of 154 spiral galaxies, in or¬ 
der to detect possible effects due to internal extinction 
by dust. The selection of the sample assures that galax¬ 
ies at different inclinations are not intrinsically different, 
so that the observed correlations represent the real be¬ 
haviour of the parameter considered. All the parameters 
are derived from a bi-dimensional fitting of the galaxy im¬ 
age. We find that extinction, though small at near infrared 
wavelengths, is sufficient to produce observable effects. In 
particular the observed increase of the average disk scale- 
length and the reddening of the disk I — H colour at high 
inclinations are clear signatures of the presence of dust. 
The total Td-band disk luminosity depends little on incli¬ 
nation; on the other hand significant corrections to the 
face-on aspect are derived for the id-band central disk 
brightness and the disk scalelength. The bulge parameters 
exhibit little or no dependence on inclination. Simulations 
carried out with a simple model for an internally-extincted 
galaxy show that these results imply a central id-band op¬ 
tical depth between 0.3 and 0.5. 
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1. Introduction 


The amount of internal extinction in spiral disks has been 
the subject of a long series of paper s throughout the last 
few years, since Disney et al. (198£, DDP hereafter) and 
Valentijn ( 199d| ), questioning the widespread assumption 
of transparency of galaxies, suggested that disks might 
actually be optically thicker than previously thought. In 
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particular their claim was in disagreement with the con¬ 
clusions by Holmberg dating back to 1958| and 1975. Sub¬ 
sequent papers (e.g. Pelet ier fe Willner |l992j ; Jansen et 
al. 1994 Giovanelli et. al 1994 - G94 hereafter) tend to 
support a scenario in which galaxies are transparent in 
the outer regions and moderately opaque at the center. 
Though this is likely to be the case, the actual amount of 
extinction in spiral disks is still a matter of debate. On 
the other hand, estimating the internal extinction in disks 
and how it affects the observed light distribution is neces¬ 
sary in order to correct the observed quantities to either 
the face-on aspect or the dust-free case. Such corrections 
are fundamental in any study concerning galaxies both as 
single objects (structure and evolution), and as tracers of 
the large scale structure of the universe. 

A possible approach to characterizing the extinction 
properties of disks relies on the statistical properties of 
large samples of galaxies, for example, analyzing the cor¬ 
relations of observable parameters with inclination. In this 
way, average corrections to the face-on aspect for such 
quantities can be calculated in principle without intro¬ 
ducing any modeling of the extinction processes. On the 
other hand, simulations of extincted galaxies (e.g. Byun 
et al. 1994 Bianchi et al. 199(j ; Corradi et al. 1996 ) can 
predict the trends observed, and therefore provide insights 
into the parameters of the dust layer, such as the thick¬ 
ness or the central optical depth r(0). A few caveats need 
to be considered when undertaking a study of this kind. 
First of all, simulations have shown that some tests are 
not significant at any optical depth, since in some cases 
a large variation in r(0) has only a minor influence on 
a particular parameter. Using many different correlations 
at the same time will necessarily yield a more confident 
estimate of the opacity. Of course every structural param¬ 
eter (scalelength, surface brightness, the inclination itself) 
needs to be determined in a most reliable and objective 
way, possibly separating the two contributions from bulge 
and disk to the total brightness distribution. If simulations 
are used to predict the observed trends with inclination, a 
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comparison is possible only if real and simulated data are 
analyzed in a similar way. In fact different decomposition 
techniques often yield discrepant estimates of the mea¬ 
sured para meter s as pointed out, e.g., by Knapen & van 
der Kruit (1992). Finally, in the galaxy sample used for 
the analysis, possible observational biases need to be ac¬ 
counted for, as demonstrated by several authors (Burstein 
et al. 1991; Davies et al. 1993; Giovanelli et al. 1995| - G95 
hereafter). Davies et al., for example, showed that galaxies 
from a diameter-limited sample tend to be characterized 
by a preferred value of the central disk surface brightness 
/i(0), which is d eterm ined by the visibility function of that 
sample (Davies 1990| ). This preferred value does not de¬ 
pend on the inclination or on the intrinsic distribution of 
surface brightness, so that any conclusion about the disk 
opacity derived by simply investigating the observed cor¬ 
relation of /i(0) with inclination is in general affected by 
this selection effect. 

In this paper, we study how the bulge and disk struc¬ 
tural parameters correlate with inclination in the near in¬ 
frared (NIR), in order to detect and quantify the effect 
of internal extinction at these wavelengths. To this pur¬ 
pose we use a sample of 154 spiral galaxies in the Pisces- 
Perseus supercluster region, for which H band images are 
available. These data were collected for a project aimed to 
investigate the structural, photometric and dynamic prop¬ 
erties in the NIR of a large homogeneous sample of spirals, 
for which a reliable estimate of the amplitude of internal 
extinction is necessary. This study can also provide infor¬ 
mation on the general issue of disk opacity, beyond its 
NIR perspective. In fact, if disks are optically thick in the 
B band, we expect them to be only partly opaque in the 
NIR, and simulations clearly show that some of the corre¬ 
lations with inclination (in particular the correlations of 
central brightness and scalelengths) are useful and signif¬ 
icant only in the case of moderate opacity. 

The separation of the bulge and disk contributions to 
the total light, and the subsequent determination of the 
various structural parameters of the two components, is 
carried out using the entire galaxy image, instead of a 
brightness profile, which considerably impro ves the relia- 
bil ity of the results (Byun & Freeman 1995 ; Moriondo et 
al. 1998 ). Since the sample is selected from a single super¬ 
cluster the scatter introduced by distance selection biases 
are minimized. For about 70 galaxies of the sample I band 
brightness profiles are also available, from the Scl sample 
of Haynes et al. (in preparation), for these objects I — H 
colour profiles were also used to investigate the extinction 
properties of the disks. 

We first describe the data analysis and discuss the cor¬ 
relations between the various parameters and the inclina¬ 
tion, summarizing briefly the trends expected at different 
opacity levels. Next we account for possible selection ef¬ 
fects on the sample and define an intrinsically homoge¬ 
neous set of galaxies. We finally discuss the trends ob¬ 
served in our sample and we compare our results with 


the predictions of a model for an extincted galaxy, de¬ 
riving an estimate for the central optical depth in the H 
band. Throughout the paper we assume a Hubble constant 
Hq = 75 km s” 1 Mpc -1 . 


2. The Sample 


Our galaxies belong to a sample selected in the region of 
the Pisces-Perseus supercluster and were initially drawn 
from the target of the redshift survey of this region car¬ 
ried out by R.G. and M.H. Selection criteria, observations, 
data reduction and photometric calibration are described 
in Moriondo et al. (in preparation). One hundred seventy 
eight galaxies were imaged in the NIR photometric bands 
(J, H , and K ), with types ranging from SO to Sd and 
covering the whole range of inclinations - though with 
fewer objects at high inclination. Since for a large frac¬ 
tion of these objects only H band data are available, in 
the remainder of the paper we will consider only the H 
band subsample (174 galaxies as a whole, mostly late-type 
spirals). 

Distances were evaluated correcting the heliocentric 
radial velocities to the CMB frame of reference with the 


prescription reported in RC3 (de Vaucouleurs et al. |l99l|) . 
The average distance of the sample is 70 Mpc with a dis¬ 
persion of 24, yielding an average FF-band absolute mag¬ 
nitude of —23.4; the luminosity range is between ~ —21 
and ~ —25 FF-mag. 


3. Data analysis 


A surface brightness profile was extracted fox each galaxy 
using the ELLIPSE routine in the STSDASB package for 
data reduction and analysis. In particular, the ellipse fit¬ 
ting was performed keeping the center of the ellipses fixed, 
but allowing ellipticity and position angle to vary. From 
such profiles we determined the isophotal radii considered 
in the following analysis. They were also used to obtain 


I — H colour profiles, as described in Sect. 6.5 


To obtain a description as reliable and objective as 
possible of the galaxy structure, rather than merely de¬ 
termining the outer slope of the major axis profiles, we 
fitted a bi-dimensional model brightness distribution con¬ 
volved with a gaussian PSF to each galaxy image. The 
model includes two components, a disk and a bulge, whose 
apparent brightness distributions are both elliptical and 
decrease exponentially with radius. In a magnitude scale, 
along the major axis: 


H(r) = At(0) + 1.086— . (1) 

rd 

1 STSDAS is distributed by the Space Telescope Science In¬ 
stitute, which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under NASA contract 
NAS 5-26555. 
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We expect an exponential shape for the bulge to be appro¬ 
priate for these galaxies since almost all of them are Sb’s 
or later types (see for example Andredakis et al. 1995 ), 
and in fact this choice is supported a posteriori by the 
good quality achieved by most of the fits. When the ex¬ 
ponential distribution for the bulge is too steep to yield 
a satisfactory fit to the data (7 cases out of 174) we con¬ 
sider a generalized distribution for this component (Sersic 
19681), i e.: 


Mb (r) = Me + 1-086 < —a„ 


1/n 

-I -i 

r, 


( 2 ) 


The exponent index n - equal to 1 for the exponential dis¬ 
tribution - is allowed to be an integer greater than 1 ; fi e 
and r e are the effective surface brightness and radius, and 
a n is a constant relating the effective brightness a nd ra - 
dius to the exponential values (see Moriondo et al. 1998 ). 
In the cases for which exponential bulge fits are unsuc¬ 
cessful a good fit to the data is obtained with n = 2. The 
parameters of each fit are the two scalelengths, the two 
surface brightnesses, and the two apparent ellipticities of 
bulge and disk. We use the fitted disk ellipticities to deter¬ 
mine the axial ratios a/b used for the correlations in the 
next sections. For 8 galaxies with particularly disturbed 
morphology we do not manage to obtain a satisfactory fit; 
these objects will therefore be excluded from the following 
analysis. More details about the fitting procedure can be 
found in Moriondo et al. 1998. 


4. Basic relations between disk parameters and in¬ 
clination 

Both the central surface brightness of the disk, /z(0), and 
its scalelength r dl are expected to correlate with inclina¬ 
tion in a way which depends on the amount of extinction 
present in the disk. This is also true for isophotal radii and 
total luminosities, which can be expressed as functions of 
fi( 0) and r d . From Eq. (!l|) a generic isophotal radius rj SO 
is defined by: 

rise = 0.921 {fiiso - /r(0 ))r d ; (3) 

whereas the total magnitude of the disk is given by 

m = ii(0) — 2.5 log(2 7 rr 2 ) + 2.5 log(a/ 6 ) . (4) 

If the galaxy is completely transparent /i(r) is expected to 
become brighter with inclination: 

A i ( r ) = /i°W - 2.5 log(a/ 6 ), (5) 

where we have denoted the face-on value with a in 
this case Eq. [l] implies that the disk scalelength will be 
independent of a/b , and from Eq. (Q) it follows that m is 
also independent of inclination, whereas from Eq. (i we 
obtain: 

T'iso Tiso 2.306 r d log(a/ 6 ) ( 6 ) 


If some extinction is present in the disk Eq. (J5j) can be 
generalized as 

M0)=ju°(0)-2.5Clog(a/&) , (7) 

where 0 < C < 1, with (7 = 0 corresponding to the com¬ 
pletely opaque case. Following G94, relationships of other 
photometric parameters with inclination are also param¬ 
eterized as linear dependences on log(a/ 6 ). For the disk 
scalelength, the isophotal radius and the total disk mag¬ 
nitude we will have respectively: 


Td 

Tiso 

m d 


= r° d [l + 77 log(a/&)] 

= r° so [t + S iso log(a/ 6 )] 
= ™° d + 7 log(a/&) . 


( 8 ) 

(9) 

( 10 ) 


Disk scalelengths are expected to be independent of in¬ 
clination (rj = 0 ) both in the case of completely transpar¬ 
ent and completely opaque disks while they will increase 
with inclination if some extinction is present in the inner 
part of the disks (Byun et al. 1994; Bianchi et al. 1996| ). 
As opacity increases, any isophotal radius will become less 
sensitive to inclination than in the transparent case (5i SO 
will decrease), whereas we expect the total luminosity to 
become a steeper function of inclination (7 will increase). 
We will assume that relations analogous to Eqs. 8-10 are 
appropriate to represent the trends of the corresponding 
bulge parameters as well. 

The relations considered so far hold both for appar¬ 
ent and absolute quantities if we consider a single galaxy. 
Yet, since we are looking for average correlations in an 
extended sample, we will use only distance-independent 
quantities, i.e. the absolute ones, for which we will main¬ 
tain a notation analogous to that shown above. 


5. Checking for selection effects 


To check for possible intrinsic differences between high and 
low inclination galaxies in our sample we assume that at 
H band galaxy disks are transparent at all inclinations at 
sufficiently large isophotal radii. In particular, we choose 
the radius where the surface brightness equals 21 H mag 
arcsec^ 2 , which corresponds on average to a distance of 
three disk scalelengths from the center and to the 22.5 
mag isophote in the I band. We note that even with a B 
band face-on optical depth as high as 20 , the triplex model 
by DDP and a standard extinction curve (e.g. Cardelli et 
al. 1989) predict an optical depth as low as 0.1 in the 
H band at about three disk scalelengths from the cen¬ 
ter. More refined models and simulations (e.g. Byun et 
al. 1994) confirm that the variation of the outer isopho¬ 


tal radii with inclination is largely independent of r( 0 ). 
In the hypothesis of transparency, we can estimate the 
face-on value of r 2 i from Eq. (||) for every galaxy in the 
sample. A fair sample should cover the same range in r%\ 
at all inclinations. 
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The quantity r° d appearing in Eq. d) is the face-on disk 
scalelength measured in the region of the galaxy where 
the disk is actually transparent. This in general will be 
different from r d , the scalelength we obtain from our de¬ 
compositions, since the latter is computed over the whole 
disk; we assume that at least for the galaxies which are 
face-on our estimate of r d can be considered a reasonable 
approximation for r° d . We correct r d for the other galaxies 
using Eq. ( 8 ), with a coefficient 77 = 0.5. This choice for 77 
for our H band data will be justified later on. The plot we 
obtain is shown in Fig. [lj here and in the rest of the paper 
Sa galaxies and earlier are represented by crosses, Sab and 
Sb by open circles, Sbc and later types by filled triangles. 
No obvious trend of r^ with inclination appears. To make 
our sample as fair as possible, spanning approximately the 
same interval in at all inclinations, we exclude 12 more 
galaxies and restrict the following analysis to the objects 
whose face-on isophotal radius lies between 2.9 and 16.9 
kpc (dotted lines in Fig. ]l[). The final sample includes 154 
galaxies. 

If we consider now the correlation of r 2 i (uncorrected 
to face-on aspect) versus log(a/ 6 ) for this subsample, we 
can obtain an estimate of the average 821 coefficient; we 
find 82 1 = 0.85 ± 0.16. From Eqs. (||) and (9) we see that 
for transparent disks, this quantity is 2.306 7 " 2 /r 21 , imply¬ 
ing r 2 ! = 2.7 r° d . On the other hand, from the observed 
ratios rd/r 2 i in the face-on case we can obtain an inde¬ 
pendent estimate for 821 '- the intercept at zero inclination 
of the best fit to r d /r 2 1 versus log(a/ 6 ) yields again 0.85 
(±0.03). Since this latter value depends also on the disk 
scalelengths derived from the image decompositions, we 
deduce that for face-on galaxies our estimate of r d is not 
significantly different from the value expected in the case 
of transparency, as we have previously assumed. 

For 151 galaxies of the final sample R.G. and M.H. 
provided 21 cm line widths, from their private database 
(“AGC”) Further evidence that our sample is not heavily 
affected by selection biases is supported by the fact that we 
do not observe any trend with log(a/ 6 ) of the line widths, 
corrected for inclination. 


6. Observed correlations with inclination 


When looking at the correlation between an observed pa¬ 
rameter and the inclination, ideally one should analyze 
separately galaxies of different morphological type and ab¬ 
solute luminosity. Objects belonging to different classes, 
in fact, can be characterized by different extinction prop¬ 
erties and therefore by different slopes in the correlation 
considered (G95, Tully et al 1998). Moreover, the observed 
parameter itself (for example the central disk brightness) 
is likely to depend both on the galaxy luminosity and the 
morphological type (G95; de Jong 1996; Morton & Haynes 


1994), and not accounting for these trends would intro¬ 


duce additional scatter in the correlation with log(a/ 6 ). In 
any correlation with inclination considered for our data, 


we do not observe obviously different slopes between galax¬ 
ies belonging to different classes (we will return to this 
point in Sect. 6.3). Therefore we assume a common slope 
for each correlation with log(a/ 6 ). On the other hand, 
as described in the next section, we account for possi¬ 
ble trends with luminosity and morphological type of the 
observed parameters, i.e. we allow the offsets of the cor¬ 
relations with inclination to be different for objects of dif¬ 
ferent classes. Such trends are rather evident, at least in 
some cases: more luminous galaxies are characterized by 
larger r d s; Sa’s are more compact on average (smaller r d 
and brighter 0 )) than later type galaxies of the same 
luminosity, and so on. The iF-band luminosity, in turn, 
appe ars to depend little on inclination, as we will show in 
Sect. 3T, and in the next two sections we will consider the 
observed luminosities, uncorrected to face-on aspect. 


6 . 1 . Disk scalelength 

Figure ||i shows the correlation between iJ-band disk 
scalelengths (on a logarithmic scale) and absolute disk lu¬ 
minosity. In order to account at the same time for the de¬ 
pendence of r d both on luminosity and on morphological 
type, the scalelengths are plotted after scaling each mor¬ 
phological type by an additive term, so that for every type, 
the average Mh and the average rescaled log(rd) lie on the 
best linear fit to the whole (rescaled) sample. The plot is 
obtained iteratively: first we perform a linear regression to 
the observed values, after rebinning the points in groups 
of 12 ; then, for each morphological type, we scale the r d s 
by an additive term so that the average scalelength for 
that class lies on the linear regression. The best fit is then 
recomputed and a new rescaling accomplished; the process 
is repeated until convergence. The residuals with respect 
to the final best fit are plotted versus log(a/ 6 ) in Fig. [l|b. 
The average scalelength is added to each residual, in order 
to preserve the same y-axis scale; the best fit of Eq. 8 to 
the plot is also shown, yielding a slope 77 = 0.41 ± 0.17. 
We remind that a correlation of r d with inclination is to 
be found only if some extinction is present, at least in the 
inner regions of the disk. 

When all the galaxies are fitted, instead of only those 
in the selected “fair sample”, we find a slope of 0.5. This 
is the value we used in Sect. || to correct the isophotal 
radii to their face-on value. We find that for any value of 
77 from 0.0 (no correlation of r d with inclination) to 0.6 
(the value derived by G94 in the I band) the plot of r 21 
versus log(a/ 6 ) (Fig. [)]) and the subsequent selection of a 
“fair sample” do not change much. As a consequence, if we 
plot the disk scalelengths versus inclination for a sample 
selected assuming 77 = 0 , the value for 77 that we derive a 
posteriori from the plot is not zero but 0.56, implying that 
an increase of r d with inclination is the only self-consistent 
possibility. 
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Fig. 1. Face-on corrected isophotal radius at 21 H mag arcsec 2 versus the logarithm of inclination. The galaxies selected for 
the following analysis are the ones between the two horizontal dotted lines (respectively at 2.9 and 16.5 kpc). Crosses correspond 
to Sa galaxies and earlier types, circles to Sab and Sb types, triangles to Sbc’s and later types. 


Fig. 2. a: Disk scalelength plotted versus absolute luminosity, after rescaling log(rd) for each morphology class; the straight 
line is the best fit to the data, b: residuals from the best fit in panel a plotted versus log(a/fo); the average scalelength is added 
to each residual in order to preserve the same y-axis scale. The best fit to r<i versus log(a/fe) is also shown, c and d: the same 
correlations for the central disk surface brightness. 


6.2. Disk surface brightness 

In principle, the behaviour of the central disk surface 
brightness /i( 0 ) is rather sensitive to the amount of ex¬ 
tinction. In practice, some authors (e.g. G94; Byun et al. 


1994) have questioned the actual reliability of the observed 


trends involving /i( 0 ) for several reasons: the large intrin¬ 
sic scatter, difficulties arising in the evaluation of /i( 0 ) due 
to the presence of the bulge, a shallow dependence of /r( 0 ) 
on inclination unless optical depth is very low. In our case, 
the use of a 2 d decomposition of the brightness distribu¬ 
tion and the fact that in the NIR extinction is actually 
lower than at optical bands allow us to attribute more 
significance to this test. 

Indeed, we do find a correlation between 0) and 
log(a/ 6 ), after accounting for possible trends of the disk 
central brightness with respect to luminosity and morpho¬ 
logical type in the way explained in the previous section. 
The plots are shown in Fig. [2|c and d. The best fit to the 
data in panel d yields C = 0.60 ± 0.14, a rather high in¬ 
crease with inclination which tends to exclude high disk 
opacities. This value is consistent with two published es¬ 
timates of the coefficient, both derived using average disk 
surface brightnesses from aperture photometry. In partic¬ 
ular Peletier & Willner ( 1992 ) confined Ch between 0.6 
and 0 . 8 , and deduced a central face-on optical depth of 
about 0.2 using a DDP sandwich model. Boselli & Gavazzi 
( 1994 ) found Ch = 0.65, independent of morphological 
type. The face-on central surface brightness turns out to 
be 17.6 H -mag arcsec -2 with a dispersion of 0.7. 


6.3. The Holmberg test 

Let us consider now the so-called Holmberg test (Holm¬ 
berg 1958| , 1975), i.e. the correlation between inclination 
and the average surface brightness within the circular 
aperture defined by a certain isophotal radius. An ob¬ 
served trend of this quantity with inclination could be 


interpreted as the signature of significant opacity in the 
disks. However, as stated by G94, a trend is expected also 
in the case of transparent galaxies because the isophotal 
radius is an increasing function of inclination (except in 
the case of completely opaque disks), so that the surface 
brightness tends to decrease at high inclinations because 
it is estimated using a larger aperture. 

Figure js| (top panel) displays the correlation found for 
our data; in this case, each morphology class was scaled 
by an additive term such as to have a common average 
on the y axis. Also shown are the expected slopes in the 
case of a completely transparent disk (dotted arrow), a 
completely opaque one (solid arrow), and the behaviour 
expected when isophotal radius and disk scalelengths vary 
with inclination according to the relations described in 
the previous sections. The three slopes are computed for a 
pure exponential and assuming r^/r^ = 2.7. We note that 
the coefficients derived from our data yield a trend which 
is nearly the same, in terms of average surface brightness, 
to the one expected for completely transparent galaxies. 
Actually, as already stated by G94, this kind of correla¬ 
tion is not particularly sensitive to the amount of internal 
extinction, and the scatter in the plot is too large to dis¬ 
tinguish reliably between the different cases. 

Following G94, we also consider a modified version of 
the Holmberg test in which the surface brightness is cal¬ 
culated within the face-on value of the isophotal radius. 
In particular, in the bottom panel of Fig. we plot the 
disk magnitude plus 5 log r%i versus the inclination, and 
the best fit to the data. Since r%\ ~ computed using Eq. 
9 with ^21 = 0.85 - is by definition independent of incli¬ 
nation, the slope of this plot is the coefficient 7 in Eq. 
( 10 ); the value of 7 corresponds roughly to the correc¬ 
tion for the total magnitude from edge-on to face-on as¬ 
pect. We find little correlation between the two quanti¬ 
ties (7 = 0.16 ± 0 . 22 ), which, at first sight, is the signa¬ 
ture of nearly complete transparency. Actually, this small 
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Fig. 3. Top panel: Average surface brightness within ?' 2 i versus log(a/6). The three arrows represent the trend expected for 
different opacities: the solid arrow is for completely opaque galaxies, the dotted arrow for completely transparent ones, the 
dashed arrow is derived from the observed slopes of I'd and fi( 0) versus inclination. Bottom panel: average surface brightness 
within r^i versus log(a/6), with the best fit to the data points. 


correlation is consistent with the other trends observed: 
the brightening of /z( 0 ) with inclination in our galaxies is 
smaller than in the transparent case, yet the increase of 
Td ~ which does not occur if disks are optically thin - con¬ 
tributes to raise the value of the total luminosity within 
r 2 i, making the average surface brightness behave very 
much like in the case of complete transparency. Quanti¬ 
tatively, if we parametrize the correlation between rd and 
(a/ 6 ) as a power law, instead of using Eq. 8 , in particular 
if we write 



it can be easily demonstrated that 7 = 1 — C — 2a. For 
our data we find a = 0.2 ± 0.05, yielding 7 = 0 ± 0.15, 
consistent with the value obtaine d from Fig. [| 

In recent work, Tully et al. ( 199S ) derived estimates 
for the extinction correction to face-on both in the optical 
and in the NIR [K' band), as a function of galaxy lumi¬ 
nosity. Assuming H — K’ ~ 0.2 mag, and since the median 
absolute luminosity for the galaxies in our sample is —23.7 
H— mag, we derive from their Eq. ( 6 ) 7 k’ = 0.26, a value 
which is consistent with ours, if we assume the difference 
between 7 k> and 7 h to be negligible. A dependence of the 
amount of internal extinction on the galaxy luminosity is 
also found by G95. For our galaxies we do not observe 
any trend of this kind, but we note that our sample cov¬ 
ers a limited range in absolute magnitude (about three 
magnitudes, and even less when considering the galaxies 
at high inclination). Moreover, for galaxies as luminous as 
ours (we estimate our /-band absolute magnitudes to be 
brighter than —20), the G95 data suggest only a rather 
shallow dependence of extinction on luminosity (see e.g. 
their Fig. 7). 


6 . 4 ■ Bulge parameters 

We look now for correlations of the //-band bulge param¬ 
eters with inclination, in particular of the bulge effective 
radius r e , effective surface brightness He and absolute mag¬ 
nitude Mb, using the same approach outlined for the disk 
parameters in the previous sections. We find a small in¬ 
crease of r e with inclination, and little or no correlation 
for fi e and Mb- In particular, if we define three coefficients 
rib, Cb and 7 & using the same relations adopted for the disk 
(Eqs. 8-10), we derive from the best fits to the data the 


following values: 

r/b = 0.24 ±0.16 (12) 

C b = 0.04 ±0.16 (13) 

7 b = 0.09 ±0.55 . (14) 

The bulge-to-disk ratio as well exhibits no significant de¬ 
pendence on inclination. 

In this case a straightforward interpretation of the re¬ 
sults is more difficult than for the disk parameters: in fact 
the bulge and the dust are characterized by completely 
different spatial distributions, respectively a peaked and 
compact spheroid and an extended thin disk, so that 
in general extinction will not be uniformely distributed 
throughout the spheroid. In the transparent case, how¬ 
ever, we would expect the scalelength to be independent 
of inclination, and the effective surface brightness to in¬ 
crease with log(a/ 6 ) (though not as much as /r( 0 ), due to 
the lower intrinsic ellipticity of the bulge). Since this is 
not the case for our sample, we interpret the different ob¬ 
served behaviour, again, as the effect of some extinction, 
that needs to be present at these wavelengths, at least in 
the central galaxy regions. 

6.5. I — H colour profiles 

For 68 of our galaxies, / band brightness profiles are avail¬ 
able from the Scl sample of Haynes et al. (in preparation); 
for these objects, we compute / — H colour profiles. To re¬ 
duce them to a common radial scale, we normalize the 
distance from the center with the value of the disk scale- 
length, measured in the H band and corrected to face-on 
aspect using Eq. ( 8 ). 

A reddening of the observed colours with inclination 
at a given radius and for a fixed stellar population is ex¬ 
pected if (and only if) some dust is present, since the opti¬ 
cal depth along the line of sight increases with inclination, 
and at shorter wavelengths the scattering is more effective 
in removing photons traveling within the plane of the disk. 
Since the dust is concentrated towards the center of the 
galaxy it will also produce a radial colour gradient, mak¬ 
ing the outer regions of the galaxy look bluer. However, 
a contribution to such a gradient can also be provided 
by a variation of metallicity or average stellar population 
with radius. Similarly, higher extinction, metallicity, and 
different stellar content, can be invoked to explain the 
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Fig. 4 .I — H radial colour profiles averaged in three bins of absolute magnitude. Radii are normalized with the disk scalelength 
corrected to face-on aspect. 


redder colours of high luminosity galaxies (Gavazzi, 1993| ; 
Tully et al. 1998). On the other hand, a correlation of 
colour with inclination provides direct and unambiguous 
evidence for reddening by dust. Figure [| shows the aver¬ 
age I — H radial profiles for our sample in three bins of 
absolute magnitude. As expected, they become bluer with 
increasing radius, and brighter galaxies tend to be redder 
at any given radius. 


We inspect for trends of I — H with inclination in three 
radial bins: within one disk scalelength, between one and 
two, and between two and three r d s respectively. We ex¬ 
clude the region inside 1/3 of ry, where different seeing 
between / and H and the presence of a bulge can produce 
extremely high ’’local” colour gradients. Since the redden¬ 
ing of colour with increasing luminosity would introduce 
scatter in the correlation with inclination, in each radial 
bin we obtain a best fit to the colour-luminosity relation 
and plot the residuals versus the inclination. We do not 
observe any difference in colour between early and late- 
type spirals; therefore in this case we didn’t rescale the 
different morphological types. Figure || (left hand pan¬ 
els) shows the observed colour-luminosity relations in the 
3 radial bins with a best fit for each bin; the slopes are 
—0.12 ± 0.01, —0.11 ± 0.02 and —0.08 ± 0.03 respectively; 
the sign of the slopes is the one expected (redder colour 
with increasing luminosity). In the right hand panels we 
plot the residuals of these correlations versus log(a/ 6 ). In 
each panel the average colour for that particular bin is 
added to each residual, to preserve the same y-axis scale; 
the best fits to the data are also shown. As expected we 
find that more inclined galaxies tend to be redder. The 
slopes of the best fits are 0.34 ± 0.04 mag, 0.56 ± 0.07, and 
0.6 ± 0.1 respectively. 


If we consider the average disk colours, the slope in 
the correlation between colour and inclination amounts to 
7 / — Jh, and we find this slope to be 0.46 ±0.14. Again, we 
can compare our result to the values found by Tully et al. 
(1998) and G94. The median I — H for our sample is 1.74, 
yielding for the I band luminosity a value of about — 22 . 
Equation (5) in Tully et al. then gives 7 / = 1.05, which 
is also the value found by G94 for their sample. Our data 
seem to suggest a slightly lower value for 7 j (i.e. lower 
extinction), for which we can set an upper limit of about 
0.9 from the colour-inclination slope and from the value 
of 7 h derived in the previous section. 


6 .6. Color gradients and inclination 

We also measured colour gradients by fitting to each pro¬ 
file a straight line, still excluding the points inside 0.3 r° d . 
Figure |(] shows the gradients obtained for this subsample 
plotted versus inclination. Assuming that an exponential 
disk fits well both the I and H brightness distributions, 
it can be easily shown that 


d(I-H) 

dp 


1.086 



(15) 


where we have denoted with ry and r# the disk scale- 
lengths in the two bands, and p = r/rn is the radius in 
units of the H band disk scalelength. If we parametrize the 
dependence of r d on inclination with Eq. [ll], we obtain 


d{I — H) 
dp 



(16) 


The fact that a negative gradient is almost always de¬ 
tected implies that on average rn is smaller than ry; more¬ 
over, this difference is about the same at all inclinations - 
amounting to 10 ~ 20 % - which is what we expect if the 
disk scalelengths correlate in a similar way with inclina¬ 
tion in both bands, i.e. an — olj. If any correlation is to 
be found in Fig. || it is in the sense of smaller gradients for 
high inclinations, implying an > a/ from Eq. ©• Look¬ 
ing at Fig. 7 in G95 this seems actually to be the case 
for galaxies brighter than —20 I- mag, for which a value of 
ay ~ 0.1 is appropriate, whereas we found an = 0.2. Note 
that from the same figure it turns out that the dependence 
on inclination of disk scalelengths is steeper (a is higher) 
for galaxies of low luminosity, which the authors claim to 
be less extincted. The fact that a is higher also in the 
H band where extinction is certainly lower is therefore a 
confirmation of G95’s result. A best fit to the data points 
in Fig. |] (dashed line) yields an — a/ = 0.11 ± 0.04, and a 
rather accurate estimate for the ratio r° H /r 0 j = 0.81 ±0.03. 


7. Disk opacity in the NIR 

Some qualitative statements can be made on the basis of 
the correlations shown above. It appears that some of the 
parameters which describe the brightness distribution of 
a galaxy are affected by the presence of interstellar dust 
also in the NIR. At these wavelengths, however, the ef¬ 
fect of extinction on most observable quantities is quite 
small. The best evidence for the presence of extinction is 
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Fig. 5. Left hand panels: colour-luminosity relations. In each panel, from top to bottom, we plot I — H colours averaged 
in a different radial bin. The best fit to the data are also shown. Right hand panels: the residuals from the best fit to the 
colour-luminosity relations in the left hand panels are plotted versus the inclination. In each panel the average colour for that 
particular bin is added to each residual, to preserve the same y-axis scale. The best fits to the data are also shown. 


Fig. 6. Radial colour gradients versus log(a/b). The dashed line is the best fit to the data assuming a power-law dependence of 
rd on the axial ratio a/b both in I and H. 


in the correlation between disk scalelength and inclination 
- which turns out to be simil ar t o the one existing in the 
I band (Sect. |6.6| and Sect. S.l ) and in that between 
I — H colour and inclination. On the other hand, the be¬ 
haviour of the disk central brightness is not very different 
from what we expect in the case of transparency; such 
a correlation, as pointed out by Byun at al, 1994 , is to 
be seen only when opacity is low. Moreover, we find little 
variation of the total luminosity with inclination, and the 
correction from edge-on to face-on for the total magnitude 
is about 0.15 mag on average confirming that extinction 
in the NIR , though detectable, is not conspicuous even 
at high inclinations. 


7.1. Simulating absorption by dust 

In order to better quantify our estimate of optical depth 
in the disks, we simulate a set of brightness distributions 
affected by absorption by dust for a range of inclinations 
and different values of the central optical depth. We adopt 
a model galaxy in which dust and stars are distributed 
exponentially both in the radial and “vertical” direction, 
like in the triplex model by DDP, but we assume the dust 
s calele ngth to be larger than rd, following Peletier et al. 
( 1995 ) , the results of two recent works by Xilouris et al. 
( 1997 , 1995 ), and of Nelson et al. 1998 . In particular our 
dust disks are 1.5 times larger than the stellar ones. The 
models are major-axis brightness profiles computed ac¬ 
cording to the approximate analytical solution for radia¬ 
tive transfer in a dusty disk described in DDP. We con¬ 
sider pure disks (no bulge component) with dust-to-stars 
scaleheight ratio ( = 0.25, 0.5, or 0.75, and r#(0) < 2. 
Scattering effects are not included. The profiles are fit¬ 
ted with exponentials excluding the inner part, within 0.3 
disk scalelengths, i.e. the region in our galaxies where the 
bulge is usually the brighter component. To make the fit¬ 
ting process as similar as possible to the one we performed 
on real galaxies, a proper weighting is given to the various 


points of the profile in order to simulate the typical errors 
in our data. Since our sample is not uniformly distributed 
in log(a/6), we define a set of inclinations for the mod¬ 
els by rebinning the values for the real data in groups of 
12 galaxies each. The highest inclination bin is centered 

O 

at 80 , where the numerical approximations in the model 
are still reasonably accurate (see DDP and G94). 

We also compute a set of models to simulate the bulge 
brightness profiles in the region where they are usually 
brighter than the disk. For our galaxies the bulge is in 
general well represented by an exponential, and it domi¬ 
nates the brightness distribution at any inclination within 
approximately two effective radii (i.e. about three bulge 
scalelengths rb, or 0.3 disk scalelengths). The median in¬ 
trinsic ellipticity is 0.4. A satisfactory modeling of the 
bulge-dominated region of a galaxy should in principle be 
carried out using two dimensional brightness distributions, 
to analyze the information provided by the asymmetries, 
due to extinction, between the two halves of the galaxy 
image. It should also account for the presence of the un¬ 
derlying disk distribution. In our case we will consider our 
simple model of the major-axis surface brightness profile 
only as a rough check for the observed trends of the bulge 
parameters with inclination. In this perspective we adopt 
the same type of model we used for the disk (the product 
of two exponentials, along the radial and vertical direction 
respectively), though a density distribution characterized 
by ellipsoidal isophotal surfaces would be more appropri¬ 
ate to represent the spheroidal shape of the bulge. For 
this set of models we choose a dust scalelength 13 times 
larger than rb] using the typical intrinsic bulge ellipticity 
(0.4) and the bulge-to-disk scalelength ratio (~ 0.1) of our 
sample galaxies, we also transform the dust-to-stars scale- 
heights adopted for the disk, namely 0.25, 0.5, and 0.75, 
into 0.35, 0.7, and 1 respectively. 
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7. ,2. Disk and bulge parameters 


Figure ^ shows the expected slopes of the correlations of 
disk scalclength, central brightness, and total luminosity 
with inclination (coefficients 77 , C, and 7 ). The coefficients 
are plotted as a function of central optical depth. In each 
panel the dashed line, the solid line, and the dotted line 
correspond to the simulations with £ = 0.25, 0.5 and 0.75 
respectively. The shaded areas are the estimates of the co¬ 
efficients derived from our sample at 1 <7 confidence level. 

It turns out that all the observed coefficients are consis¬ 
tent with the models in a range of central optical depths. 
In particular, from panel b, we can set a lower limit for 
th( 0) around 0.3 regardless of the thickness of the dust 
layer, since the three coefficients considered exhibit lit¬ 
tle dependence on £ for th( 0) ^ 0.5. From panel c, on 
the other hand, we derive an upper limit to the optical 
depth around 0.5 again roughly independent of £. We note 
that our lower limit implies a larger th ( 0 ) than the values 
around 0.1 suggested by the results in Xilouris et al. (1997, 


1998). Such values, however, are derived for two particu¬ 


lar edge-on galaxies and not from a statistical analysis, 
and we do not consider the discrepancy particularly wor¬ 
risome. Adopting the extinction curve reported by Gordon 
et al. 1997 for the Milky Way, our limits for th{ 0) imply 
a central optical depth between 0.9 and 1.4 in the / band, 
between 1.8 and 3 in the V band, and between 2.4 and 
3.9 in the B band. In agreement with most of the recent 
works on the subject, therefore, our data suggest that spi¬ 
ral disks are characterized by a moderate opacity also at 
optical wavelengths. In particular the presence of a sig¬ 
nificant extinction in the outer disk regions can be safely 
ruled out, except in the case of extremely inclined galax¬ 
ies. On the other hand, the extinction corrections derived 
from our data cannot be easily extrapolated to the B or 
V band on the basis of the extincton curve alone, since 
they also depend on the exact geometry of dust and stars, 
and on the wavelength- depen dent effect of scattering (see 
for example Byun et al. 1994 ). 

In our model we assumed a dust disk larger by 50% 
than the stellar one. If we adopt equal scalelengths for 
the two components we find that the predicted slopes are 
consistent with our data for slightly higher values of th{ 0 ), 
namely between 0.4 and 0.8. If we consider instead an 
even more extended dust layer, in particular if the dust 
scalelength in twice as large as r^, again we find agreement 
between models and data for around 0.3. 


Figure H is the analogous of Fig. |i] for the bulge param¬ 
eters. Comparing the model predictions in the two figures 
the largest difference is to be found in the behaviour of the 
coefficient C, which for the disk is equal to 1 in the trans¬ 
parent case and then decreases tending to zero as r( 0 ) 
increases. For the bulge the apparent axial ratio is always 
> b/a, so that in the transparent case Cb < 1 ; moreover, 
since the dust distribution is much more extended than 
the bulge in the radial direction, for sufficiently high opti¬ 


cal depths the surface brightness along the major axis be¬ 
comes fainter with increasing inclination, yielding Cb < 0. 

As explained in the previous section these simulations 
are to be considered less reliable than the ones carried out 
for the disks, and actually there are no values of the cen¬ 
tral optical depth for which we have agreement between 
models and data for all the three correlations considered, 
at 1 a confidence level. We note however that at 2 er any 
th( 0) ^ 0.5 is consistent with the data; as was the case 
for Fig. [ 7 ] this is true regardless of the dust-to-stars scale- 
height ratio. This range of optical depths, in turn, is con¬ 
sistent with the one derived from the disk correlations, 
so that we can consider the results from these tests as a 
confirmation - though at a lower confidence level - of our 
previous findings. 


7.3. I — H colours 


Figure |] shows the dependence on th{ 0) of the simulated 
slopes of / — H colour versus log(a/6) in the three radial 
bins considered for the real data. The values we derive 
from our galaxies are represented by the shaded areas. In 
the simulations 77 is assumed to be 2.87 times the value in 
the II ba nd, according to the value reported by Gordon et 
al. ( 1997 ) for the Milky Way. From this plot, we infer that 
the slopes derived from our data are consistent with the 
limits previously set on th{ 0) (0.3 ~ 0.5) only if £ < 0.5. 
For thicker dust layers the dependence of the colour on 
inclination tends to be too steep within r^, unless we in¬ 
crease the central optical depth above 1. This result holds 
regardless of the radial extension of the dust distribution 
with respect to the stellar scalelength. A thick dust distri¬ 
bution appears unlikely also because of the little difference 
observed in the dependence of rd on inclination between 
the H and I band, i.e. for two optical depths differing by 
about a factor of three (see Sect. 6.6). From Fig. ^a, we 
see that 77 is constant or decreasing with increasing r only 
for low values of £ (we neglect differences in £ between H 
and I: e.g., Xilouris et al. 1997 find this ratio to be fairly 
constant from K to B for NGC 2048). 


7 . 4 . Is scattering important? 

In our simple model for an internally-extincted galaxy, we 
have neglected the effect of scattering of radiation by dust; 
roughly speaking, this means that for a given amount of 
dust (i.e. for a given r(0)) our model galaxy is less lumi¬ 
nous than a real galaxy with the same structural parame¬ 
ters. How does this affect the conclusions derived from our 
simulations? We note that we are not interested on how 
scattering affects the observed quantities, but on how it 
affects their correlations with inclination. A few Monte 
Carlo simulations of a dusty spiral galaxy kindly provided 
by S. Bianchi allowed us to quantify approximately the 
differences introduced by scattering in the observed 77- 
band brightness distribution at low and high inclinations. 
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Fig. 7. From top to bottom: 17-band r/, C, and 7 coefficients versus central optical depth. Dashed lines represent simulations 
with dust-to-stars scaleheight ratio £ = 0.25; for thin solid lines £ = 0.5, and for dotted lines £ = 0.75. The shaded areas are 
the ranges we derive from the observed disk parameters at 1 a confidence level. 


Fig. 8. From top to bottom: 17-band r/b, Cb, and jb coefficients versus central optical depth. Dashed lines represent simulations 
with dust-to-stars scaleheight ratio £ = 0.35; for thin solid lines £ = 0.7, and for dotted lines £ = 1. The shaded areas are the 
ranges we derive from the observed bulge parameters at 1 <7 confidence level. 


Fig. 9. Slope of the correlation of 7 — 77 colour versus inclination, plotted versus central 77-band optical depth, in three radial 
bins. Dashed lines represent simulations with dust-to-stars scaleheight ratio £ = 0.25; for thin solid lines £ = 0.5, and for dotted 
lines £ = 0.75. The shaded areas are the ranges we derive from our data at 1 a confidence level. 


The code which performed the simulations, described in 
Bianchi et al. ( 1996 ), produces a set of images at different 
inclinations for a given galaxy computing the radiative 
transfer without introducing any approximation and in¬ 
cluding both absorption and multiple scattering by dust. 
For every simulated galaxy image the corresponding image 
including only absorption was also computed. A compar¬ 
ison of the two sets of images (with and without scatter¬ 
ing) for a galaxy with the average structural properties of 
our sample and a central optical depth 0 ) = 1 reveals 
that the contribution by scattering to the observed sur¬ 
face brightness is slightly higher at low inclinations and 
increasing toward the center. Yet, the implied corrections 
to rd and /r( 0 ) are rather insensitive to inclination when 
compared to the observed variations of these parameters 
with log(a/ 6 ). We estimate the correction to the central 
disk brightness due to scattering to be roughly 0.06 77- 
mag arcsec ^ 2 larger in the face-on case than in the case 

O 

of a highly inclined galaxy (~ 80 ), and the decrease in 
rd to be around 2% larger in the face-on case. Both these 
corrections are negligible, since the observed correlations 
with inclination for these parameters imply a brighten¬ 
ing of ~ 1.15 mag for /r(0), and an increase of ~ 30% 

O 

for rd, between face-on and 80 inclination. The differ¬ 
ence in the scattering correction to the total magnitude is 
more significant (~ 0.04 mag) given the small variation of 
this quantity with log(a/ 6 ), but still rather small. More¬ 
over, since the central optical depth is likely to be less 
than 1 in the 77 band, the actual corrections are probably 


even smaller than these. The effect of scattering is likely 
to be more important when considering the correlations 
with inclinations of the bulge parameters, given the little 
dependence on log(a/ 6 ) we find for all of them, and the 
highly concentrated bulge brightness distribution. On the 
other hand, in this case our model is a rather poor approx¬ 
imation anyway, and we considered its predictions only as 
a rough confirmation of the results derived from the disk 
parameters. 

In conclusion, we expect that the inclusion of scatter¬ 
ing in the model would not produce significant changes in 
our estimate of t/j ( 0 ). 


8. Summary 

In this paper, we have derived NIR (77 band) struc¬ 
tural parameters for a sample of 174 spiral galaxies in 
the Pisces-Perseus supercluster, using a bi-dimensional de¬ 
composition of the galaxy images. We have subsequently 
analyzed the correlations of such parameters with incli¬ 
nation in order to detect possible effects of internal ex¬ 
tinction and estimate the relative corrections to face-on 
aspect. For 68 galaxies of the sample 7 — 77 colour profiles 
are also available, which were used to the same purpose. 
The main results of this work are: 

1) The effect of internal extinction in the 77 band, 
though small, can be detected. In particular the increase of 
the average disk scalelength and of the disk 7 — 77 colours 
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with inclination can be attributed to the presence of the 
dust. 

2) We find little dependence of the total disk lumi¬ 
nosity on inclination. We derive an average correction of 
about 0.15 mag from edge-on to face-on aspect in the H 
band. We find this result to be consistent with the corre¬ 
lations we find for and /r(0) versus inclination. Little 
or no correlation with inclination is found for the bulge 
paramjeters. 
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3) Comparing our results with the predictions of a sim¬ 
ple model for extincted galaxies, we deduce that the av¬ 
erage central optical depth in the H band is between 0.3 
and 0.5 if the dust scalelength is 50% larger than r^. The 
simulations also suggest that the ratio of the dust-to-stars 
scaleheights is on average lower than 0.5. The effect of 
scattering on the various correlations considered is found 
to be negligible at these wavelengths. 


Xilouris E.M., Kylafis N.D., Papamastorakis J., et al., 1997, 
A&A 325, 135 

Xilouris E.M., Alton P.B., Davies J.I., et al., 1998, A&A 331, 
894 

Valentijn E.A. 1990, Nat 346, 153 
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